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I ! Airborne  satellite  communication  systems  have  been  shown  to  frequently 
suffer  severe  degradation  of  performance  due  to  ionospheric  scintillation. 
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The  objective  of  this  study  was  to  lay  the  groundwork  for  a quantitative 
evaluation  of  the  performance  improvement  that  can  be  achieved  through  the 
use  of  an  adaptive  error  control  technique.  Specifically,  the  adaptive  Galiager 
algorithm  and  the  extended  Galiager  algorithm  were  examined  to  determine 
which  of  these  provides  the  better  solution  to  the  error  control  problem  on 
this  channel. 

The  first  phase  of  this  study  dealt  with  the  heuristic  analysis  of  these  two 
error  control  techniques  on  the  basis  of  binary  error  sequences  generated 
from  flight  test  records  of  received  signal  level.  Phase  I culminates  with 
a recommendation  as  to  which  technique  should  be  used.  Suggestions  re- 
garding the  values  of  important  algorithm  parameters  are  also  provided. 


Phase  II  of  this  study  was  devoted  to  the  development  of  a computer  simu- 
lation algorithm  and  the  design  of  a hardware  evaluator  for  the  recommended 
oodlng  technique.  The  algorithm,  used  with  the  AFAL  channel  simulator,  permits 
the  quantitative  evaluation  of  code  performance  on  the  simulated  scintillation 
channel. -The  hardware  evaluator,  which  is  an  adaptive  Galiager  encoder/ 
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switch  selectable  parameters,  is  designed  to  be  interface 
Ifh  the  existing  AFAL  flight  test  modem. 
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SECTION  I 
INTRODUCTION 

| I 

Airborne  satellite  communication  systems  have  been  shown  to  frequently 
suffer  severe  degradation  of  performance  due  to  ionospheric  scintillation.  This 
type  of  interference  produces  deep  fades  in  received  signal  power  resulting 
in  long  error  bursts  in  the  demodulated  data  stream. 

The  objective  of  this  study  was  to  lay  the  groundwork  for  a quantitative 
evaluation  of  the  performance  improvement  that  can  be  achieved  through  the 
use  of  an  adaptive  error  control  technique.  Specifically,  the  adaptive  Gallager 
algorithm  and  the  extended  Gallager  algorithm  were  examined  to  determine 
which  of  these  provides  the  better  solution  to  the  error  control  problem  on 
this  channel. 

The  first  phase  of  this  study  dealt  with  the  heuristic  analysis  of  these 
two  error  control  techniques  on  the  basis  of  binary  error  sequences  generated 
from  flight  test  records  of  received  signal  level.  Phase  I culminates  with  a 
recommendation  as  to  which  technique  should  be  used.  Suggestions  regarding 
the  values  of  important  algorithm  parameters  are  also  provided. 

Phase  II  of  this  study  was  devoted  to  the  development  of  a computer 
simulation  algorithm  and  the  design  of  a hardware  evaluator  for  the  recommended 
coding  technique.  The  algorithm,  used  with  the  AFAL  channel  simulator,  permits 
the  quantitative  evaluation  of  code  performance  on  the  simulated  scintillation 
channel.  The  hardware  evaluator,  which  is  an  adaptive  Gallager  encoder/ 
decoder  with  switch  selectable  parameters,  is  designed  to  be  interface  com- 

* 

patible  with  the  existing  AFAL  flight  test  modem. 
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SECTION  II 


PHASE  I - ANALYSIS  OF  CODES 
1.  INTRODUCTION 

The  objective  of  Phase  I of  this  study  was  to  determine  heuristically, 
based  on  error  records  supplied  by  AFAL,  whether  the  adaptive  Gallager  or 
the  extended  Gallager  error  correcting  algorithm  offers  the  best  solution  to 
the  error  control  problem  for  this  ionospheric  scintillation  channel.  In 
doing  this,  a search  was  conducted  for  new  convolutional  codes  with  the 
properties  necessary  for  use  with  the  extended  algorithm.  This  search  was 
unsuccessful  in  that  no  new  high  rate  convolutional  codes  with  the  properties 
required  for  use  with  the  extended  algorithm  were  discovered.  (The  extended 
algorithm  requires  two  convolutional  codes;  one  containing  the  other.) 

The  search  for  new  convolutional  codes  was  conducted  using  the 
trial -and -error  technique  with  the  aid  of  a threshold  decoding  simulator  that 
was  developed  for  this  purpose.  Although  no  new  "contained"  codes  were 
discovered,  two  trial  polynomials  were  confirmed  as  generators  for  1/2  rate 
convolutional  codes  with  good  distance  properties. 

2 . ADAPTIVE  GALLAGER  ALGORITHM 

A diagram  of  a typical  1/2— rate  Gallager  encoder  is  shown  in  Figure 
1.  It  consists  of  a shift  register  B+X+k  bits  long  with  several  taps  at 
the  left  end  and  one  tap  on  the  right-most  stage.  The  configuration  of 
taps  at  the  left  end  of  the  register  is  defined  by  the  code  generator  poly- 
nomial which  is  chosen  based  on  good  random  error  correction  and  detection 
properties. 
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Figure  1 . Adaptive  Gallager  Encoder/Decoder 


As  each  information  bit  is  shifted  into  the  register,  it  is  transmitted 
over  the  channel.  In  the  next  channel  signalling  interval,  a parity  bit  is 
generated  and  sent  by  forming  the  modulo  2 sum  of  the  oldest  bit  i,  and 
the  new  bits  in  the  tapped  stages  at  the  left  end  of  the  register. 

The  adaptive  Gallager  decoder  is  illustrated  in  Figure  lb.  The 
decoder  contains  a replica  of  the  encoder  shift  register.  The  parity  generator, 
however,  has  as  an  additional  input  the  received  parity  bit.  Thus,  the 
modulo  2 adder  generates  a parity  check  on  the  received  information  bits  and 
compares  it  to  the  corresponding  received  parity  bit  to  produce  a syndrome 
bit.  Examination  of  the  syndrome  sequence  gives  information  about  the  error 
pattern. 

The  decoder  has  two  modes  of  operation  between  which  it  switches 
automatically  based  on  the  channel  condition.  If  the  channel  condition  is 
such  that  the  random  error  correcting  code  can  handle  the  errors , error  correction 
is  performed  at  stage  ix  of  the  decoder  register  as  directed  by  the  Random  Error 

Corrector  logic.  The  decoder  criterion  operates  on  the  syndrome  sequence; 

(1) 

typically  a threshold  decoder  of  the  Massey  type  is  used.  If  the  random 
error  correcting  power  of  the  code  is  not  exceeded,  errors  will  be  corrected 
before  reaching  stage  ij  . The  decoder  tap  at  stage  ij  will  thus  have  no 
effect  on  the  syndrome  since  the  data  sequence  is  transparent  to  the  syndrome. 

The  principle  of  operation  of  the  decoder  is  based  on  the  ability  of 
the  code  to  detect  certain  error  patterns  of  greater  weight  than  it  can  correct. 
When  an  uncorrectable,  but  detectable,  error  pattern  occurs,  the  random 
error  correction  control  algorithm  is  suspended  and  burst  error  correcting 


takes  over.  The  burst  correcting  algorithm  is  very  simple.  The  rule  is  to 
decide  that  ij  is  wrong  and  change  it  if  and  only  if  the  most  recent  syndrome 
bit  SD.  ...  =1.  If  the  error  burst  is  no  more  than  2B  bits  long,  it  willhave 

passed  the  taps  at  the  left  end  of  the  decoder  register  before  the  decoder 
enters  the  burst  mode.  If  a "clean"  guard  space  follows  the  error  burst,  the 
current  syndrome  will  be  non-zero  (flagging  a correction)  if  and  only  if  ij 
itself  is  in  error.  If  a guard  space  of  error-free  bits  approximately  equal  to 
the  burst  length  is  present  at  the  input  to  the  decoder,  the  entire  error  burst 
will  be  decoded  by  this  rule.  Although  not  shown  in  the  diagram,  when  a 
correction  is  made  the  effect  of  the  error  on  the  syndrome  is  removed  by 
complimenting  SB+x+)c  (making  it  a binary  zero).  Because  of  this,  as  the 
error  burst  passes  out  of  the  decoder,  the  syndrome  register  fills  up  with 
zeros.  The  control  algorithm  examines  a region  around  the  right  end  of  the 
syndrome  register  for  an  all  zero  condition.  When  enough  consecutive  zero 
syndromes  are  observed,  control  is  passed  back  to  the  random  error  correcting 
algorithm. 

The  advantage  of  the  adaptive  Gallager  algorithm  over  many  of  the 
other  burst  correcting  codes  is  that  it  requires  a very  short  error-free  guard 
space  approximately  equal  in  length  to  the  actual  burst  that  is  being  corrected. 

Other  techniques  typically  require  a guard  space  about  three  times  the  length 
of  the  maximum  correctable  burst.  The  reason  for  this  is  that  the  Gallager 
technique  corrects  most  but  not  100%,  of  the  bursts  less  than  its  maximum 
designed  length.  The  requirement  for  a guard-space-to-burst  ratio  of  3 applies 
to  the  idealized  case  of  guaranteed  error-free  burst  correction. 
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Note  that  the  occurrence  of  an  error  in  the  guard  space  will  cause  a 
short  burst  of  errors  in  the  output  (one  each  time  it  passes  through  a tapped 
stage  of  the  register)  if  the  error  is  in  an  information  bit  and  a single  output 
error  if  the  error  is  in  a parity  bit.  This  problem  can  be  more  or  less  serious 
depending  on  the  error  statistics  on  the  channel  of  interest.  For  the  scintil- 
lation channel,  the  random  error  rate  in  the  interburst  intervals  is  rather  high. 

For  this  reason,  the  effect  of  guard  space  errors  is  an  important  consideration. 

One  approach  to  dealing  with  the  problem  of  guard  space  errors  is  to 
modify  the  rules  that  govern  the  burst/random  mode  switching.  Indications  as 
to  the  modifications  that  may  be  required  can  be  obtained  by  running  computer 
simulations  using  the  standard  algorithm  and  an  error  sequence  generator 
whose  statistics  accurately  model  the  channel  of  interest.  The  simulation 
should  log  the  occurrence  of  decoded  bit  errors  and  keep  track  of  the  state 
of  the  decoder  (burst  or  random)  when  each  of  these  failures  occurs.  Analysis 
of  the  results  of  such  a simulation  will  give  clues  as  to  the  failure  modes  of 
the  algorithm.  For  example,  if  a preponderance  of  the  errors  occur  in  the 
burst  mode,  the  criteria  for  entering  the  burst  mode  may  need  to  be  strengthened 
so  that  the  decoder  doesn't  switch  to  the  burst  mode  so  readily.  On  the  other  hand, 
weakening  the  requirements  for  switching  back  to  the  random  mode  may  be  the  answer 
The  encoder  and  decoder  simulation  algorithms  provided  in  Section  III  of  this 
report  can  be  used  in  conjunction  with  the  AFAL  channel  simulation  program 
to  conduct  an  analysis  such  as  this.  The  hardware  evaluator  design  presented 
in  Section  III  provides  for  a limited  amount  of  experimentation  in  this  regard 
through  the  provision  of  a switch  that  allows  the  random-to-burst  criterion  to  be 
strengthened.  In  addition,  the  design  permits  the  use  of  any  desired  algorithm 
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(stored  in  programmable  read-only  memory)  to  control  the  basic  random-to- 
burst  switching  decision.  Computer  simulation  can  be  used  to  determine  the 
optimum  ROM  program  for  the  scintillation  channel.  Finally,  the  ease  with 
which  the  decoder  returns  to  the  random  mode  can  be  affected  by  varying  the 
'Y*  - parameter. 

3 . EXTENDED  GALLAGER  ALGORITHM 

While  the  performance  of  the  adaptive  Gallager  decoder  can  be  optimized 

for  specific  channels  by  employing  the  proper  mode-change  strategies,  there  is 

a definite  limit  to  the  performance  that  can  be  achieved  in  this  manner.  Another 

approach  to  the  problem  of  errors  in  the  guard  space  has  been  suggested  by 

(2) 

Sullivan.  Sullivan's  generalized  (or  extended)  Gallager  decoder  depends , in 
its  principle  of  operation,  upon  the  use  of  two  convolutional  codes,  C and  C*, 
where  C*  contains  C.  At  the  encoder  the  information  sequence  is  first  encoded 
using  code  C;  after  a fixed  delay,  it  is  also  encoded  with  a "shortened"  version 
of  C*,  which  is  added  to  the  parity  bits  of  C . A functional  diagram  of  the 
extended  Gallager  decoder  is  shown  in  Figure  2.  In  the  random  error  correcting 
mode,  the  decoder  is  equivalent  in  operation  to  the  ordinary  Gallager  decoder 
using  the  C-  code.  In  the  burst  mode,  the  C*  decoder  is  switched  in  to  re- 
move random  errors  from  the  interval  following  the  error  burst. 

In  his  paper,  Sullivan  presented  an  example  of  the  extended  Gallager 
algorithm  using  convolutional  codes  that  yield  an  equivalent  coding  rate  of 
1/3.  He  also  points  out  the  difficulty  in  the  general  application  of  his 
scheme  by  warning  that  there  is  "a  fundamental  problem  . . . resulting 


from  the  scarcity  of  good  constructive  random-error-correcting  convolutional  codes, 
particularly  at  high  rates.  This  problem  becomes  even  more  pronounced  with 

the  constraint  that  one  of  the  two  codes  used  must  contain  the  other."  He  goes 
on  to  say,  "It  therefore  appears  likely  that  full  utilization  of  this  scheme 
must  await  further  developments  in  constructive  techniques  for  the  encoding 
and  decoding  of  random-error-correcting  convolutional  codes." 

Subsequent  to  the  appearance  of  Sullivan's  paper,  several  papers 
which  bear  on  the  problem  have  appeared.  In  May,  1972,  a correspondence 
entitled,  "Contained  Convolutional  Codes,"  appeared  in  the  Transactions 
on  Information  Theory  (Ferguson, 3)  . Here,  Ferguson  shows  that  the  conditions 
of  containment  imply  a very  simple  factorization  of  the  codes.  More  recently, 

Wu (4)  , (5)  in  a two-part  paper,  has  presented  a code-generation  algorithm  for 
high  rate  convolutional  codes.  Since  these  codes  are  threshold  decodable, 

Wu's  algorithm  has  potential  application  to  the  problem  of  discovering  good 
codes  for  use  with  the  Gallager  algorithm.  Unfortunately,  the  appearance 
of  Wu's  paper  was  subsequent  to  the  completion  of  Phase  I of  this  study. 

4.  SEARCH  FOR  NEW  CODES 

During  the  conduct  of  Phase  I of  this  study,  M.  Kim  of  ECI  conducted 
a search  for  new  convolutional  codes  with  the  properties  required  for  ap- 
plication to  the  extended  Gallager  algorithm.  This  search  was  carried  out 

(1) 

using  the  trial  and  error  method  together  with  a threshold  decoding  simulation 
that  was  developed  to  test  the  properties  of  codes  generated  by  the  various 
trial  generator  polynomials.  A listing  of  this  computer  simulation  is  presented 
in  Appendix  A. 
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The  impetus  for  this  attempt  to  find  a pair  of  more  powerful  convo- 
lutional codes  for  use  with  the  extended  Gallager  algorithm  came  from  the  fact 
that  the  known  codes  of  rate  1/3  are  not  powerful  enough  to  correct  3 consecutive 
errors  in  the  random  mode  or  to  correct  fc/vo  consecutive  errors  in  the  guard 
space  in  the  burst  mode. 

Although  the  search  for  contained  codes  was  unsuccessful,  Kim  did 
succeed  in  discovering  two  new  polynomials  that  generate  codes  with  good 
distance  properties.  These  are  shown  in  Table  1 where  the  properties  of 
some  of  the  known  codes  as  well  as  the  new  codes  are  summarized.  The  new 
code  of  constraint  length  24  was  selected  for  the  proposed  hardware  evaluator 
design. 

5.  CHANNEL  ERROR  CHARACTERISTICS 

A communication  channel  can,  in  general,  be  categorized  into  one  of 
three  basic  classes:  random  error  channels,  burst  error  channels  and  compound 
channels  (combination  random  and  burst).  The  UHF  channel  affected  by 
ionospheric  scintillation  fading  falls  into  the  compound  channel  category. 

Eight  files  of  measured  received  signal  level  taken  from  actual  flight 
tests  together  with  binary  error  sequences  generated  by  a computer  model  of 
the  modem  were  furnished  to  ECI  to  provide  a basis  for  an  heuristic  analysis 
of  code  performance. 
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TABLE  1 . CANDIDATE  CONVOLUTIONAL  CODES 


FOR  ORIGINAL  GALLAGER  DECODING  (1/2-RATE) 

GENERATING 

EFFECTIVE 

GUARANTEED 

POLYNOMIAL 

CONSTRAINT 

CORRECTABLE 

DETECTABLE 

LENGTH  (N£) 

ERRORS 

ERRORS 

REMARKS 

(111101) 

16 

2/20 

3/20 

(111001) 

11 

2/12 

3/12 

(100000110111) 

22 

3/24 

3/24 

(101101110111) 

24 

3/24 

4/24 

New  Code 

(10110111) 

16 

2/18 

3/18 

New  Code 

FOR  EXTENDED  GALLAGER  DECODING  (1/3 

-RATE) 

(2) 

Gj  =(1101) 

g/3)  = (1100) 

— 

2/12 

3/12 

C-Code 

G2*(3)=  (1100) 
*(3) 

Gx  = (0111) 

« __ 

1/12 

2/12 

C*-Code 
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The  following  is  a brief  synopsis  of  the  AFAL  program  that  was  used 
to  analyze  the  channel  error  statistics. 

The  flight  data  file  number  and  position  on  disk  are  read  in  from  the 
keyboard.  To  approximate  a 75  BPS  transmission  rate,  a variable  is  set 
so  as  to  use  only  every  14th  data  point  from  flight  file.  This  is  arrived  at 
by  knowing  that  each  file  is  approximately  5 minutes  long.  There  being  301 
blocks  of  data  in  a file  makes  a single  block  correspond  to  one  second.  Each 
block  contains  1066  data  points  which  roughly  corresponds  'o  14  X 75  (1050). 
The  data  points  are  not  averaged  (14  points  at  a time)  because  the  digitized 
power  level  was  very  smooth.  A printout  of  a digitized  block  is  included  in 
the  material.  Seven  burst  threshold  values  were  investigated  which  are  .15, 
.12,  .10,  .08,  .05,  .03,  .01.  Each  burst  threshold  value  was  evaluated 
at  4 different  .001  threshold  values  that  correspond  to  -136  dbm,  -134  dbm, 
-132  dbm,  -130  dbm.  The  minimum  burst  length  is  set  at  1/4  second.  Burst 
and  guard  counts  cleared  and  flags  initialized. 

The  program  now  evaluates  burst  length  versus  guard  length  for  300 
blocks  of  a flight  data  file.  The  following  criteria  are  used  for  burst  and 
guard  analysis. 

The  detection  of  a burst  is  anytime  the  bit  error  probability  exceeds  or 
equals  the  burst  threshold  level.  This  generally  signifies  the  end  of  a 
guard  space.  All  burst  counts  below  1/4  second  are  padded  to  at  least 
a 1/4  second  burst  by  part  of  the  guard  space  to  the  burst.  All  burst 
counts  equaling  4 seconds  duration  are  counted  and  the  burst  count  is  put 
back  to  one.  This  is  done  because  a guard  count  equaling  1/4  of  the  burst 
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length  may  not  arise  so  that  useful  information  v/ould  be  lost.  A burst  count 
is  complete  whenever  the  guard  count  is  at  least  equal  to  1/4  of  the  burst 
length.  The  burst  lengths  for  each  file  are  tabulated  in  1/4  second  increments 
.25  - .50,  .50  - .75,  etc.  up  to  3 seconds  with  burst  between  3 and  4 seconds 
in  one  group  and  those  over  4 seconds  counted  separately.  A guard  count  is 
started  whenever  the  bit  error  probability  falls  below  the  burst  threshold  level 
and  the  previous  burst  length  being  at  least  1/4  second.  A guard  count 
terminates  whenever  a burst  condition  occurs  (anytime  bit  error  probability 
exceeds  or  equals  burst  threshold  level) . Whenever  a guard  count  is  less  than 
1/4  of  the  burst  count,  the  guard  count  is  added  to  the  burst  count.  The  guard 
lengths  for  each  file  are  tabulated  corresponding  to  1/4  of  the  previous  burst 
length  .25  - .50  BL,  .50  - .75  BL,  etc.  up  to  3 times  the  burst  length.  All  guard 
lengths  greater  than  3 times  length  are  tabulated  together.  When  all  300  blocks 
of  data  (corresponding  to  22,500  transmitted  bits)  have  been  processed,  the 
tabulated  data  is  printed  out  on  the  lineprinter  and  written  out  on  the  disc. 

Among  the  8 files  of  received  signal  level  that  were  provided  by  AFAL, 

File  No.  4 appears  to  represent  the  worst  case.  Tabulated  burst  statistics  for 
both  the  aggregate  of  all  8 files  and  File  No.  4 are  shown  in  Table  2.  The 
conclusions  presented  in  the  following  section  are  based  primarily  on  analysis 
of  the  binary  error  sequences  generated  from  File  No.  4 data. 

5.  CONCLUSIONS 

Analysis  of  binary  error  sequences  generated  from  File  No.  4 data 
reveals  that  the  longest  error  burst  is  639  bits  in  length.  Most  of  the  long 
error  bursts  are  separated  by  guard  spaces  less  than  2B  bits  in  length.  The 
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TABLE  2 . BURST  ERROR  STATISTICS 
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total  number  of  error  burst  occurrences  in  the  21 ,500  bit  record  was  27.  The 
error  burst  duty  cycle  was  approximately  33%  with  an  average  burst  length 
of  263  bits.  While  the  average  guard  space-to-burst  ratio  was  approximately  3, 
for  long  bursts,  the  G/B  ratio  was  close  to  unity.  The  background  error  rate 
was  moderately  low.  However,  there  was  a relatively  high  incidence  of 
short  error  clusters  consisting  of  2 or  3 consecutive  errors. 

Because  of  the  high  incidence  of  consecutive  errors  in  the  guard  spaces, 
it  is  concluded  that  the  performance  improvement  that  would  result  from  the 
use  of  the  extended  algori  thm  would  not  be  significantly  better  than  the 
benefit  that  can  be  derived  from  the  original  Gallager  algorithm.  This  is  because 
the  known  extended  Gallager  algorithm  codes  of  rate  1/3  are  not  powerful  enough 
to  correct  3 consecutive  errors  in  the  random  mode  or  2 consecutive  errors  in 
the  guard  space  in  the  burst  mode.  Furthermore,  the  extended  algorithm 
requires  1.5  times  more  transmission  time  and  nearly  twice  the  hardware 
complexity  of  the  original  Gallager  algorithm.  Based  on  this  heuristic  reason- 
ing, it  was  concluded  that  Phase  II  of  this  program  should  concentrate  on 
the  development  of  software  and  hardware  tools  for  the  quantitative  evaluation 
of  the  performance  of  the  ordinary  adaptive  Gallager  algorithm  on  the  UHF 
ionospheric  scintillation  channel.  The  new  convolutional  code  of  constraint 
length  24  specified  by  the  generator  polynomial  given  in  Table  1 should 
provide  a significant  improvement  in  decoded  error  rate  for  this  channel.  Since 
this  code  is  capable  of  correcting  3 consecutive  errors  in  the  random  mode,  this 
power  coupled  with  a burst/random  mode  change  strategy  that  is  optimized 
for  the  scintillation  channel  should  approach  the  performance  that  could  be 
achieved  with  the  more  costly  1/3  rate  extended  Gallager  algorithm. 


15 


of  course,  adjustable  in  both  the  simulation  algorithm  and  the  hardware 
evaluator. 


SECTION  III 


PHASE  II  - DESIGN  OF  EVALUATION  TOOLS 

1 . INTRODUCTION 

As  a result  of  Phase  I of  this  program,  it  was  determined  that  the 
expected  performance  improvement  Drovided  by  the  extended  Gallager 
algorithm  on  the  scintillation  channel  does  not  warrant  the  substantial 
added  complexity.  Because  of  this.  Phase  II  of  this  program  has  been 
devoted  to  the  development  of  tools  for  evaluating  the  performance  of  the 
original  Gallager  algorithm  on  a real  scintillation  channel.  Specifically, 
software  algorithms  tailored  to  the  PDP-11  have  been  developed  for  use 
with  the  AFAL  channel  simulation  program  for  the  adaptive  Gallager  encoder 
and  decoder.  In  addition,  a detailed  hardware  design  has  been  completed 
for  a variable  parameter  adaptive  Gallager  encoder-decoder  that  will  enable 
real-channel  evaluation  of  the  effectiveness  of  this  error  correcting  scheme. 

Description  of  these  software  and  hardware  performance  evaluation 
tools  is  the  subject  of  this  section. 

2.  GALLAGER  ALGORITHM  USING  g(D)  = 1+D2+D3 *+D5+D6+D7 *+D9+Dl0+Dn 

Block  diagrams  of  the  coder  and  decoder  for  the  code  selected  are 

shown  in  Figures  3 and  4 respectively.  The  encoder  is  extremely  simple, 

consisting  of  a single  long  buffer  register,  a parity  tree  and  a commutator 

that  alternately  sends  data  and  parity  bits.  As  designated  by  the  generator 

polynomial,  parity  bits  are  computed  from  certain  of  the  past  twelve  in- 

formation bits  together  with  the  information  bit  in  the  last  stage  of  the  buffer. 

Since  this  is  a 1/2  rate  code,  two  channel  symbols  (one  information  and  one 
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Figure  3.  Block  Diagram  - Adaptive  Gallager 
Encoder  (Kim's  Polynomial) 


parity)  are  sent  for  each  information  symbol  that  enters  the  encoder. 


The  decoder  contains  a replica  of  the  encoder  register  and  parity 
tree  plus  the  additional  computational  and  decision  logic  needed  to  im- 
plement the  decoding  algorithm.  This  consists  in  a second  B+X+ll  stage 
syndrome  buffer  and  two  decision  algorithms  that  operate  on  certain  segments 
of  the  syndrome  register.  (The  syndrome  register  stores  the  sequence  of 
syndrome  bits.  Each  syndrome  bit  is  generated  by  comparing  the  received 
parity  bit  with  a regenerated  parity  bit  computed  from  the  received  information 
bits).  Based  on  the  content  of  stages  X through  X+ll,  one  of  these  algorithms 
decides  in  favor  of  one  of  three  alternatives  each  time  the  register  is  ad- 
vanced. These  alternatives  are,  "Do  nothing,"  "Correct  the  bit  in  stage 
'X'  of  the  decoder  register",  or  "Switch  the  decoder  to  the  burst  mode  of 
operation."  On  channels  which  do  not  exhibit  clean  guard  space  between 
error  bursts,  it  is  desirable  to  prevent  entering  the  burst  mode  if  errors 
are  present  at  the  decoder  input.  In  order  to  implement  this,  logic  can  be 
added  to  prevent  entering  the  burst  mode  if  there  is  evidence  of  errors 
in  the  recent  past.  The  dashed  line  in  the  block  diagram  indicates  this 
option  which  results  in  conditioning  the  third  alternative  of  the  above 
algorithm  to,  "Switch  the  decoder  to  the  burst  mode  provided  the  most 
recent  'N'  bits  of  the  syndrome  are  all  'zero'." 

The  second  algorithm  operates  on  the  last  'Y*  stages  of  the  syndrome 
register.  This  is  a simple  rule  that  results  in  either  the  instruction,  "Do 
nothing"  (if  the  last  'Y‘  stages  contain  any  "ones")  or  the  instruction, 

"Switch  the  decoder  to  the  random-error  mode"  (if  "all-zeros"  are  present 
in  the  last  'Y'  stages). 
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3.  ALGORITHMS  FOR  PDP-11  SIMULATION 

Because  of  the  fact  that  the  AFAL  channel  simulation  was  implemented 
using  encoder  and  decoder  subroutines  written  in  PDP-11  assembly  language, 
the  algorithms  provided  here  are  structured  to  simplify  their  implementation 
using  PDP-11  registers  and  instructions.  The  schematic  drawing  shown  in 
Figure  5.  depicts  the  representation  of  the  encoder  and  decoder  buffers  in 
the  PDP-11  memory. 

Although  the  encoder  and  decoder  algorithms  are  specified  here  in 
the  form  of  considerably  detailed  How  charts,  the  tasks  of  coding  and  inter- 
facing these  subroutines  with  the  AFAL  simulation  program  remain.  (The 
coding  for  the  major  portion  of  the  encoder  subroutine  is  included  here  as 
an  example  of  one  way  the  long  buffer  registers  can  be  implemented  and 
manipulated  within  the  framework  of  the  PDP-11.) 

Flow  charts  of  the  encoding  and  decoding  algorithms  are  shown  in 
Figures  6 and  7 respectively.  These  algorithms  are  general  specifi- 
cations of  the  encoding  and  decoding  rules.  In  the  interest  of  simplifying 
the  implementation  on  the  PDP-11,  the  programmer  may  find  it  convenient 
to  impose  some  restrictions  on  the  lengths  of  the  'B'  and  'X'  registers. 

This  has  been  done  to  some  extent  in  the  example  encoder  program  given 
in  Table  3.  Here  the  encoder  register  is  represented  in  PDP-11  memory 
by  a block  of  contiguous  words.  This  approach  (rather  than  using  one 
memory  word  per  bit)  restricts  the  buffer  length  to  a multiple  of  16.  This  is 
of  no  great  consequence  performance-wise,  however,  on  the  scintillation 
channel  where  bursts  of  errors  are  typically  much  longer  than  16-bits. 
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CARRY  FLIP-FLOP 


Figure  5.  Representation  of  Decoder  Registers 
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Figure  6,  Flow  Chart  - Adaptive 
Gallager  Encoder  Simulation 
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Figure  7.  Flow  Chart  - Adaptive 
Gallager  Decoder  Simulation 
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TABLE  3 . EXAMPLE  ENCODER  PROGRAM 
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TABLE  3.  EXAMPLE  ENCODER  PROGRAM  (CONTINUED) 
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The  decoding  algorithm  presented  here  is  specialized  to  the  extent 
that  the  conventional  majority  decision  approach  of  using  a set  of  orthogonal 
equations  as  the  decision  variable  is  adopted.  If  desired,  however,  the 
algorithm  can  easily  be  modified  to  accomodate  any  alternate  decision 
rule  for  making  the  correction/mode  change  decision. 

The  decoding  algorithm  includes  the  option  of  conditioning  the  random 
to  burst  mode-change  decision  on  the  requirement  that  there  have  been  no 
errors  at  the  decoder  input  in  the  "recent  past"  as  suggested  by  Forney  ^ . 

The  purpose  of  this  is  to  prevent  entering  the  burst  mode  too  easily  since  this 
can  result  in  an  increase  in  the  decoded  bit  error  rate  on  channels  that  exhibit 
diffuse  rather  than  dense  error  bursts.  This  option  is  invoked  by  setting  the 
'DIFFUSE'  flag  equal  to  one. 


4 . HARDWARE 

The  literature  on  error  correcting  codes  is  replete  with  testimony  and 
words  of  caution  with  regard  to  the  pitfalls  and  inaccuracies  that  can  be 
encountered  in  attempts  to  predict  error  control  code  performance  on  the 
basis  of  channel  models.  It  is  almost  universally  agreed  that  the  "proof 
of  the  pudding"  can  only  come  through  real-channel  testing. 

Since  AFAL  has  the  facilities  for  scintillation  channel  flight  testing 
and  a test  modem  that  currently  contains  a feedback  encoder/decoder  card, 
it  was  proposed  that  ECI  design  a flexible  adaptive  Galiager  encoder/de- 
coder with  compatible  interfaces.  Accordingly,  a design  has  been  produced 
for  a variable-parameter,  adaptive  Galiager  decoder.  With  the  implementation 
described  here,  all  of  the  encoder  and  decoder  parameters  affecting  performance 


can  be  varied  over  a wide  range.  With  the  exception  of  the  low-power 
Schottky  read-only  memory,  the  entire  codec  has  been  designed  with  CMOS 
logic. 


5.  ADAPTIVE  GALLAGER  ENCODER. 

A schematic  diagram  of  the  encoder  is  shown  in  Figure  8.  This 
Figure  shows  that  the  Gallager  encoding  algorithm  is  simple  to  implement 
in  hardware.  Parity  bit  generation  is  implemented  with  a single  MSI  package 
(a  CMOS  12-bit  parity  tree)  while  the  parity  and  data  are  commutated  to 
form  the  output  data  stream  using  an  and-or  select  circuit.  An  output  clock 
at  twice  the  input  clock  frequency  is  generated  using  a combination  of 
exclusive-or  gates  and  a one-shot. 

A normal  design  implementation  of  the  Gallager  encoder,  including 
output  clock  generation,  would  require  only  7 I.C.’s  if  an  MOS  LSI  serial 
register  were  used  for  the  encoder  buffer.  Here,  however,  interest  lies  in 
a flexible  test  codec  with  parameter  variability.  For  this  reason,  the 
encoder  buffer  has  been  designed  with  two  variable-length  sections.  The 
primary  variable  buffer  is  labeled  R-l  in  Figure  8.  A schematic  diagram 
of  the  implementation  of  the  buffer  using  three  CMOS  LSI  Quad  64-bit  shift 
registers  (Fairchild  34731's)  is  shown  in  Figure  g,  As  shown  in  the  table 
accompanying  the  Figure,  this  implementation  provides  a very  wide  selection 
of  buffer  lengths.  The  buffer  length  selector  switch  is  ganged  to  the  cor- 
responding decoder  buffer  length  selectors. 

The  need  for  the  second  (8-stage)  variable  register  arises  from  the 
decoder  design.  In  addition  to  the  basic  buffer  length  parameter,  the  decoder 
has  two  other  design  parameters  that  affect  the  operation  of  the  device.  One 


Figure  8.  Schematic  Diagram  - Variable  Parameter  Adaptive 
Gallager  Encoder  (Kim's  Polynomial) 


Figure  9.  Variable  Buffer  Implementation 
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Figure  10.  Schematic  - Variable 
Parameter  Gallager  Decoder 


of  these  parameters  (denoted  as  the  X-parameter)  is  related  to  the  integration 
time  allotted  to  the  random-to-burst  mode  switch  decision.  The  X-parameter 
is  equal  to  one  greater  than  the  number  of  stages  of  delay  between  the 
random  error  and  burst  error  correction  points  in  the  decoder  buffer.  Since 
the  X-parameter  is  variable  and  affects  the  overall  buffer  length,  the  8-stage 
variable  section  is  included  in  the  encoder  buffer  with  the  selector  switch 
ganged  to  the  corresponding  switch  in  the  decoder  in  order  to  maintain  equal 
buffer  lengths  between  encoder  and  decoder. 


6.  DECODER  IMPLEMENTATION 

A circuit  diagram  of  the  variable  parameter  Gallager  decoder  is  shown 
in  Figure  10.  With  the  exception  of  a bipolar  PROM,  the  design  uses  CMOS 
logic  exclusively.  In  the  following  description  of  operation,  reference  is 
made  to  the  control  signals  defined  in  the  timing  diagram  shown  in  Figure  11 

The  input  serial  bit  stream  is  fed  to  one  input  of  a 12-bit  parity  tree 
in  addition  to  the  input  stage  of  the  decoder  buffer  register.  The  buffer 
register  clock,  Rj},  and  the  syndrome  register  clock,  Rgt  are  compliments  of 
each  other.  The  clock  frequency  of  each  of  these  is  one-half  the  incut  bit 
rate.  Because  of  this,  every-other-bit  of  the  input  stream  is  clocked  into 
the  data  buffer  while  the  intervening  bits  are  presented  to  the  parity  generator 
during  the  rising  edge  of  the  syndrome  register  clock.  Thus,  even  numbered 
bits  (data)  are  clocked  into  the  data  buffer  and  odd  numbered  bits  (parity) 
are  modulo  2 added  to  regenerated  parity  to  produce  a syndrome  sequence 
that  is  shifted  into  the  syndrome  register.  Upon  start-up,  the  initial 
clock  phases  may  be  such  that  parity  bits  are  clocked  into  the  data  buffer 
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INPUT  CLOCK 


Figure  11.  Timing  Diagram 


and  data  bits  are  presented  to  the  parity  checker.  This  synchronization 

problem  and  its  solution  are  discussed  in  Section  7 . 

The  following  description  of  the  sequence  of  decoder  operations  is 
illustrated  by  the  flow  chart  shown  in  Figure  12  together  with  the  schematic 
diagram  of  Figure  10  The  description  of  operation  begins  with  the  occurrence 
of  a leading  edge  of  the  syndrome  register  clock.  When  this  occurs,  the  syndrome 
register  is  right-shifted  one  bit  position.  If  the  decoder  is  in  the  BURST  mode 
and  the  current  syndrome  bit  is  a 'one',  a 'zero’  is  shifted  into  the  first  stage 
of  the  syndrome  register.  The  reason  for  this  is  explained  later.  After  a short 
delay  to  allow  for  shift  register  and  gate  delays,  strobe  STRl  samples  the  output 
of  the  NOR  gate  network  that  is  used  to  look  for  all  zeros  in  the  last  'Y'  stages 
of  the  syndrome  register.  If  an  all-zero  condition  exists,  a reset  pulse  to  the  MODE 
flip-flop  is  generated.  In  the  reset  state,  the  MODE  flip-flop  indicates  the  RANDOM 
error  correction  mode. 

The  next  event  that  occurs,  is  the  examination  of  stages  X through  X+ll 
of  the  syndrome  register  in  order  to  determine  if  a random  error  correction  is  to 
be  performed,  if  the  unit  is  to  be  switched  to  the  burst  mode  cf  operation,  or  if 
no  action  is  required.  The  rule  that  is  used  in  deciding  in  favor  of  one  of  these 
alternatives  is  defined  by  the  logic  contained  in  the  block  labeled  X-l  in  the 
schematic  diagram.  Two  alternate  implementations  of  the  block  X-l  are  shown 
in  Figure  13. 

Approach  X-la  uses  a majority  logic  decision  algorithm  based  on  a set  of 
orthogonal  equations  involving  the  syndrome  bits.  The  majority  decision  is 
performed  by  a 64  x 2 field  programmable  read-only  memory.  There  are  four 
possible  ROM  outputs,  two  of  which  are  "do  nothing"  indications.  A 
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do  nothing"  output  results  when  there  are  four  or  more  "zeros"  on  the 


address  lines.  If  there  are  four  or  more  "ones"  in  the  address,  the  ROM 
will  indicate,  "perform  a random  error  correction."  If  there  are  an  equal 
number  of  "ones"  and  "zeros"  in  the  address,  the  ROM  will  issue  the 
command  to  "switch  to  burst  mode,"  The  derivation  of  the  set  of  orthogonal 
equations  for  the  24,  12  code,  together  with  the  rationale  for  the  majority 
logic  approach  is  presented  in  Appendix  b. 

While  the  majority  logic  approach  is  the  standard  implementation  found 
in  the  literature,  the  optimum  decision  algorithm  in  any  particular  situation 
depends  greatly  on  the  channel  error  statistics.  Approach  X-lb  shows  a ROM 
implemented  decision  algorithm  where  the  (here  unspecified)  ROM  mapping 
operates  directly  on  the  syndrome  register  contents.  Extensive  computer 
simulation  using  records  of  channel  errors  would  be  required  in  order  to 
determine  the  optimum  ROM  map  for  a particular  channel. 

The  decoder  flow  chart  shows  the  three  possible  paths  taken  by  the 
decoder  as  indicated  by  the  ROM  output.  If  the  ROM  indicates  a ramdom  error 
correction  and  the  MODE  flip-flop  is  in  the  reset  state  (indicating  that  the 
decoder  is  in  the  RANDOM  mode),  the  R-CORRECT  flip-flop  will  be  set  so 
that  the  correction  will  be  effected  when  the  data  clock  edge  occurs. 

If  the  ROM  indicates  "switch  to  burst  mode,"  the  action  that  takes 
place  is  conditional,  depending  on  the  position  of  the  DIFFUSE  switch.  If 
the  DIFFUSE  switch  is  in  the  CLEAN  GUARD  position,  the  MODE  flip-flop 
is  unconditionally  set  effecting  a switch  to  BURST  mode.  If  the  DIFFUSE 
switch  is  in  the  DIFFUSE  BURST  position,  the  switch  to  BURST  mode  is 
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effected  only  if  there  have  been  no  recent  errors,  i.e. , the  last  twelve 
syndrome  bits  are  all  zero. 

If  the  ROM  gives  a "do  nothing"  output,  nothing  happens,  of  course, 
until  strobe  STR3  occurs. 

When  STR3  occurs,  the  B-CORRECT  flip-flop  will  be  set  to  the  'one' 
state  if  the  D-input  line  is  high.  This  condition  requires  that  the  decoder 
be  in  the  BURST  mode  and  that  the  current  syndrome  bit  be  a ’one’.  (The 
principle  of  operation  depends  upon  the  assumption  that  the  entire  error 
burst  is  contained  in  the  data  buffer  so  that  if  the  parity  check  fails  - indicated 
by  a syndrome  of  'one'  — this  could  only  be  caused  by  an  error  in  the  output 
stage.)  If  the  B -CORRECT  flip-flop  is  set,  the  content  of  the  output  stage 
will  be  inverted  the  next  time  the  data  buffer  is  shifted. 

The  next  event  that  takes  place  (after  the  occurrence  of  STR3)  is  the 
shifting  of  the  data  buffer  and  the  correction  of  errors  if  so  indicated  by 
the  states  of  the  correction  flip-flops . 

The  final  action  in  the  decoding  cycle  is  triggered  by  the  occurrence 
of  STR4.  The  operations  initiated  by  STR4  include  removal  of  the  effects  of 
any  errors  that  were  corrected  on  the  content  of  the  syndrome  register  and 
the  resetting  of  the  correction  flip-flops  in  preparation  for  the  next  decoding 
cycle.  In  the  case  where  a random  error  correction  has  been  performed, 
removal  of  the  effect  of  the  error  on  the  syndrome  sequence  consists  in 
complimenting  certain  bits  in  stages  X through  X+ll  of  the  syndrome  register. 
This  is  accomplished  by  incorporating  two  parallel-out  registers  that  re- 
dundantly store  syndrome  bits  Sx  through  Sx+H.  After  a random  error 
correction  is  performed,  the  contents  of  these  registers,  with  the  appropriate 

bits  inverted,  are  "jammed"  into  syndrome  register  stages  X through  X+ll. 
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The  same  general  method  is  used  to  take  out  the  error  effect  in  the 
burst  mode.  In  this  case,  however,  the  effect  is  removed  by  clearing  certain 
of  the  first  twelve  stages  of  the  syndrome  register  and  clearing  the  most 
recent  syndrome  bit.  The  former  is  accomplished  synchronous  with  STR4 
while  the  latter  is  accomplished  by  clocking  a 'zero'  into  the  syndrome  regis- 
ter whenever  a burst  correction  is  about  to  be  performed. 

The  decoder  design  provides  for  a wide  range  of  variability  of  the 
important  Gallager  decoder  parameters.  The  parameter  to  which  decoder 
performance  will  be  most  sensitive  is  the  basic  buffer  length,  B.  As  selected 
using  switches  SWl  and  SW4  (see  Figure  9)  the  buffer  length  can  be  varied 
from  12  stages  to  816  stages  in  relatively  fine-grain  increments.  The  formula 
for  buffer  length,  B,  in  terms  of  the  SW4  setting,  i,  and  the  SWl  setting,  j, 
is  given  by: 

B = 18i  + 64j  + 12  ; i = 0,3 
; J = 0,12 

Both  SWl  and  SW4  are  3-section  rotary  switches  so  that  the  encoder,  decoder, 
and  syndrome  buffers  are  varied  together. 

The  X-parameter  is  equal  to  one  greater  than  the  number  of  stages  of 
delay  between  the  points  in  the  data  buffer  at  which  random  and  burst  cor- 
rection is  performed.  The  time  the  decoder  has  in  which  to  make  a random- 
to-burst  mode  change  decision  is  thus  proportional  to  this  parameter.  Very 

little  analytical  or  empirical  information  is  available  to  guide  the  selection 

(7) 

of  this  parameter  value  (Brayer  used  a value  of  X = 15  in  his  simulations). 
Here,  the  value  has  been  given  a range  of  variation  of  from  14  to  21  under 
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control  of  selector  switch  SW2.  A simple  wiring  change  (eliminating  the 
final  4-stages  of  each  of  the  buffers)  would  shift  this  selection  to  the  range 
X = 10  to  X = 17.  Selector  switch  SW2  is  a 4 section  rotary  switch.  One 
section  is  devoted  to  each  of  the  three  buffers;  the  fourth  section  is  used 
to  disable  the  inputs  to  the  burst-to-random  mode  switch  decision  logic 
that  are  associated  with  unused  syndrome  buffer  stages. 

The  Y-parameter  is  equal  to  the  number  of  syndrome  bits  that  are 
examined  in  making  the  burst-to-random  mode  change  decision.  The  Y parameter 
sets  the  degree  of  confidence  needed  in  the  decision  that  the  burst  has  passed 
in  order  to  decide  in  favor  of  a change  to  the  random  mode.  If  Y is  made 
too  small,  errors  may  remain  in  the  data  buffer  which  will  be  passed  to 
the  output.  The  Y parameter  is  set  by  selector  switch  SW3  which  is  used 
to  disable  from  0 to  12  inputs  to  the  mode  change  decision  logic.  This 
parameter  can  thus  be  varied  from  Y = X to  Y = X+ll . 

7.  SYNCHRONIZATION 

Since  this  is  a convolutional  code,  no  real  word  sync  is  required. 

The  sync  problem  consists  only  of  determining  the  proper  phasing  of  the 
decoder  data  and  syndrome  register  clocks  so  that  data  bits,  not  parity,  are 
clocked  into  the  data  buffer. 

The  two  system  clocks  shown  in  Figure  11  are  derived  simply  by 
dividing  the  received  bit  clock  by  two  and  using  the  Q and  0 outputs  of  the 
divide-by-two  flip-flop  as  shown  in  Figure  14.  The  synchronization  problem 
arises  from  the  uncertainty  that  exists  with  regard  to  the  relationship  between 
the  incoming  bit  stream  and  the  Initial  clock  states.  If  the  start-up  state  is 


Figure  14.  Timing  and  Synchronization  Circuitry 


such  that  the  rising  edge  of  the  data  register  clock  occurs  when  a parity  bit 
is  present  at  the  input,  parity  bits  will  erroneously  be  clocked  into  the  data 
register.  This  condition  will  occur  upon  start-up  with  probability  1/2. 

Two  means  are  provided  to  cure  an  erroneous  start-up  condition.  One 
of  these,  the  manual  method,  provides  a pushbutton  switch  that  can  be 
depressed  to  invert  Doth  clock  signals  if  the  operator  observes  that  the  output 
is  'garbage.'  The  alternate  method  of  sync  acquisition  is  by  means  of  an  automatic 
sync  circuit  whose  principle  of  operation  is  based  on  the  fact  that  an  out-of-sync 
condition  will  result  in  the  occurrence  of  nearly  50%  'ones'  in  the  syndrome 
sequence . 

The  auto  sync  circuit  is  comprised  of  a monostable  with  a variable  rate 
from  300  Hz  to  6KHz,  an  'and'  gate  and  two  14-stage  ripple  counters.  The 
syndrome  sequence  is  "and-ed"  with  the  monostable  output  so  that  a series 
of  clock  pulses  is  sent  to  the  integrating  counter  each  time  a 'one'  occurs 
in  the  syndrome  sequence.  By  varying  the  monostable  frequency,  the  number 
of  pulses  produced  for  each  occurrence  of  a 'one'  can  be  varied  from  4 to  80. 

The  integrating  counter  is  periodically  reset  to  prevent  the  flagging  of  an 
"out-of-sync"  indication  during  proper  in  sync  operation.  The  frequency  of 
these  resets  is  adjustable  by  selecting  the  proper  output  tap  on  the  reset 
generating  counter.  The  auto  sync  circuit  parameters  (both  counter  outputs 
and  the  monostable  frequency)  should  be  adjusted  experimentally  so  that 
"out-of-sync"  indications  are  not  falsely  generated  when  even  the  longest 
error  bursts  are  present. 

A parts  list  for  the  hardware  evaluator  is  shown  in  Table  4. 
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TABLE  4.  ENCODER/DECODER  PARTS  LIST 


Qty . 

Part  Number 

Description 

Vendor 

1 

CD4001 

Quad  2- Input  NOR  Gate 

RCA 

2 

CD4002 

Dual  4-Input  NOR  Gate 

RCA 

6 

CD4006 

18-Stage  Shift  Register 

RCA 

1 

CD4011 

Quad  2 -Input  NAND  Gate 

RCA 

3 

CD4013 

Dual  D-Flip-Flop 

RCA 

5 

CD4015 

Dual  4-Bit  Shift  Register 

RCA 

2 

CD4020 

14-Stage  Ripple  Counter 

RCA 

7 

CD4030 

Quad  Exclusive  OR  Gate 

RCA 

4 

CD4034 

8-Stage  Bidirectional  Shift 

Register 

RCA 

1 

CD4047 

Multivibrator 

RCA 

2 

CD4049 

Hex  Inverter 

RCA 

5 

CD4078 

8- Input  NOR  Gate 

RCA 

1 

CD4082 

Dual  4-Input  AND  Gate 

RCA 

3 

CD4098 

Dual  Neonostable  Multivibrator 

RCA 

8 

CD4081 

Quad  2 -Input  AND  Gate 

RCA 

1 

CD4555 

Dual  Binary  to  1 of  4 Decoder 

RCA 

12 

74C164 

8-Bit  Parallel-Out  Shift  Register 

National 

2 

74S387 

1024-Bit  Programmable  ROM 

TI 

2 

MC14531 

12-Bit  Parity  Tree 

Motorola 

9 

34731 

Quad  64-Bit  Shift  Register 

Fairchild 

3 

851-000-U2J0- 

101J 

100  pf  Ceramic  Capacitor 

Erie 

1 

8121-100-C0G0- 

33K 

330  pf  Ceramic  Capacitor 

Erie 

1 

8121-050-651- 

103M 

.01  y f Ceramic  Capacitor 

Erie 

1 

lOOKft  1/4W  5%  Resistor 

3 

180KQ  1/4W  5%  Resistor 

1 

3099P 

lMfi  Cermet  Trimpot 

Bourns 

1 

3099P 

2 MO  Cermet  Trimpot 

Bourns 

1 

SA31SDT6 

Pushbutton  Switch 

Cutler  Hammer 

2 

SF11SCT691 

SPDT  Toggle  Switch 

Cutler  Hammer 

1 

2H50A16-3 

2-16  Position  3 Pole  Switch 

Cutler  Hammer 

1 

2E00A24-1  - 

Progressive 

Opening 

24  Position  Switch 

Cutler  Hammer 

1 

399720JC 

2-10  Position  Consecutive 

Shorting  Switch 

Oak 

1 

399475JC 

4-Section  2-12  Position  Switch 

Oak 

44 


SECTION  IV 
CONCLUSIONS 


Although  the  original  Gallager  algorithm  was  chosen  for  further  study 
on  the  basis  of  the  Phase  I analysis,  the  extended  Gallager  algorithm  should 
not  be  excluded  from  consideration  for  use  on  the  scintillation  channel  in  the 
future.  There  are  two  reasons  for  this.  First,  advances  in  coding  theory  may 
soon  yield  constructive  procedures  for  the  generation  of  more  powerful  con- 
volutional codes  with  the  containment  property  required  by  the  extended  algorithm. 
Secondly,  advances  in  integrated  circuit  technology  are  driving  the  cost  of 
digital  hardware  down  at  such  a rate  that  the  difference  in  implementation 
cor^lexity  may  soon  become  an  insignificant  consideration. 

This  study  has  produced  the  tools  necessary  for  conducting  a quantative 
determination  of  the  effectiveness  of  the  adaptive  Gallager  error  correcting 
method.  The  evaluation  can  be  approached  via  two  avenues:  computer  simu- 
lation using  the  algorithm  provided  or  flight  tests  using  the  hardware  evaluator. 
Since  the  hardware  involved  is  relatively  simple,  the  design  is  largely  com- 
pleted and  flight  tests  will  yield  the  most  meaningful  results,  this  is  the 
approach  that  is  recommended. 
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Threshold  DECODED  simulator — -m.f.  ktm  ragf  1 


//  for 

*ONF  WORD  TNTFGFRS 

SUBPOUTTNF  CO OFR . 

COMMON  L^GTH,  TLAct«NP.KR,LX.L0 
Common  TNOMT(lon) 

common  KB£fi  n = ,3)  «kgfn  US.3«&)  .MCOOF  <] 00.3) 

COMMON  MODF.  LL*  JTH,  L0con»  l°CNF 
COMMON  JONTflOO.?) . JSTRGfl 00.?) . TFPOP<300) 

COMMON  JFYVOdOO.?)  . JF  (3)  • JPOIJT  (3 ) .KOOTH (1 5 . ? . 1 S ) .KO^T!  (IF.?.  IF) 
C CONVOLUTION At  CODE  GENERATOR.  Bs  KR/NP 


T DATA  INOFX  (PIT  NUMOFR)  TL  AcT TWR  NIIMBfr  OF  DATA 

k — - SHIFT  REG  INDEX  krfG((L.k> kp  SmIFT  oEGIfTFof 

LNGTh CONFTP A I NT  PLOCK  LENGTH  j CODF  WORD  TNpfx 

mr  OUTPUT  PIT  MUmpfp  (N)  NCOOf(J.N)  — CODE  WORD  fTog 

kO IMPMTT  PIT  MUMPFR  (L) 

KGFN(N.L.K.)  CODE  OFNFPATORF 

OFFTNITION  of  FIIMC  MOOTU(M) 

MOOTtH  M ) =M_ ( M/p ) *? 

INPUT  AND  OUTPUT  POINTING  FORMAT  

4 FoRMATtlOX.  7011) 

_6_ FORMAT  (IX////  30X  . XCONVOIJITTONAI  FMCOOFR*/// 

HOX.jSOATA  TNDFXrf.  1 OX  , *f  NOUT  OATAX  . I OX  . *FNCOOFO  OATAF// 

PSOX.FGl  (I)  *.3X*0?<TM.3X,*03(T)*.3X.*G4  <!)*//) 

INITIALIZATION  

LDLAY=LB*L* 

LSYND=LX*L0*l  NOTH 
DO  1 OF  L = l.KO 
DO  1 OF  < = 1 .LNGTM 
KRFG<f.L)=0 
I OF  CONTINUE 


INC'JT  INFORMATION 

PF  AO  ( P * 4 ) ( INPUT  (T)  »I  = l.Tt_AFT) 

CONVOLUTIONAL  EMCOOTNG 

1*1 

JMR= ( TLAST/kr) 

QO  4RQ  J=t . JMR  

DO  290  L=1  .k-p 
KsLNGTH 
211  *<•  1 =K-1 

IF(K1)?20.??0.21C 
21F  KPEG (K*L) =KPEG («1 .L ) 

K«K-1 

GO  TO  211 

220  yrfG(k.l)=TNPUT (T) 

1=1*1 

?90  CONTINUE 

PARITY  CHECK  ° I Tc  

N=  1_ 

DO  390  N=1 .NR 
IF  (n-KR)  300.300.310 
310  NTEMPsO 
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THRESHOLD  OFCOnc-ts  SIMULATOR 
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LSI 

K=1 

DO  3S0  L=1»KR 
no  340  K=1*LNGTH 

MTEmp=<GENCY«L»N)  *kREG(K,L)  *MTFMP 
340  CONTINUE 
3S0  CONTINUE 

NCODFf J.N) sMODTUINTEMO) 

GO  TO  300 

300  NCODE ( J* n ) iWWFB (1 ♦ M ) 

300  CONTINUE 
400  CONTINUE 

return 

END 

//  OliR 

*<TORE  WS  UA  COOfO 

//  rno 

•ONE  WORD  INTEGERS 

subroutine  jonow 

COMMON  LNGTH,TLAOT,NR,KR,LX,LR 

COMMON  INPUT (100)  

" COMMON  KRE0(]C,3) «kgEm ( 1 S « 3 • 4 ) •NCOOE ( 1 o 6 * 3 ) 

COMMON  MOO^.  LL  • JTH « LOCOO.  (RCHC 

COMMON  JOUT< 100 *2)  » JSTRGdOO.R)  .IEPOR(300) 

COMMON  JSYWOUOO.2)  .JE<3)  .JROMT(3)  .KORTH  (15«?» 15)  « KOn T 1 (1  5 , ? , \ c > 
C INPUT  VAR  I ARLES 

C _ LOCOO-  PARITY  CHECY  length  __  _ 

C js  = NUMBER  nr  parity  eg. 

C JTH  = THREOHOLO  LEVEL 

C JOCHS  sDImmenston  OE  Rarity  EO. 

“C  KORThTK  » NS » JS ) = LINEAR  COMRINATTON  OE  RAPTTY  EO. 
moDTU(M)=m_(m/2)*P 
NS-MP-KR 
ml=LL 
JPADO=0 

no  2100  je=i*jpchs 

JORTH=0" 

DO  2000  N=1«MS 
K=LB+LMGTh 

TT-LOCOO 

2020  JOPTHsJSYNn(K»M) *KORTH ( K 1 « N « jo ) ♦JORTH 
IE(K1-1)2030*2030.?02E 
2025  K=K-1 

K1=K1-1 
GO  TO  2020 
2030  GO  TO  2000 
2000  CONTINUE 

jPAon=MonTtM  jortw) *jp Ann 
2100  "CONTINUE 

IE  (JDADn-jTH)  R110« 2120. 2130 
2110  JE (Ml  ) =0  _ 
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2120  JF(ML)*0 
JP0FT=1 

MnOE=l  

GO  TO  2200 
2130  JF(ML)=1 

J1=LR 

JTFMP= JS'TOG  ( J 1 . 1 ) ♦ 1 
jstrg(  ji  ♦md=mootim  jtfmp) 

2200  pF~UPM  _ 

end 

//  mip 

• STOOP  ws  1 1 A JPOOM 

7778B 

*ONF  WORO  imTFGPos 

• • X HRESHOLDlNG  OECODEP*  WRITTEN  o Y M.  k"TM 

• IOCS  (FARO.  1 1 3?  PRINTER, TYPEWPTTER.k'FYROARD) 

DIMENSION  IRC'/O  (100,3)  * JPGEN  (4  ) 

COMMON  LMOTH, TLAST,mo,kP.L  X,LP 

_ COMMON  fwoilT  (TOO  > 

COMMON  KRFO  ( IS  *3)  ♦ k’GEN  < 1S.3,4)  , MCOOE ( 1 00 » 3 > 

COMMON^ MOOF*  LL*_JTH»  LOCOO,  JPCHS  

COMMON  JOIJT  (100,2)  , JSTRG(100»2)  ,TFROR(300) 

COMMON  .ISYNO  ( 1 00,2)  » Jc ( 3)  » JROIJT  ( 3 ) »KOPTH ( 1 5 , 2, 1 5 ) ,KO=Tl  (IS,  3,1=:) 

OFFIN I T ION  OF  FUNCTION 

MOOTO(M)=M.(M/2)*2 
MOOUL  (M)  sM.  (M/NR)  *NR 

C INPUT  ANO  OUTPUT  FORMAT  

S FORMAT (fOX .13) 

Q FORMAT ( 1X///J  OX*  rfORTHOOOMAL  LTNPAR  COMPINATTON*//) 

J0_  F O R NAT  (10X,*P*.?Tl,*(0)=*,15n/> 

14  FORMAT  (T0*  .70T1> 

15  FORMAT (IX///// 

l iqx, /threshold  ofcootmg  simulation  for  gallagfp  cor»lN«p// 

21 X » FT  IMF  TNOEX*,3X  , *OAT A TNPUT* * 3X • /ENCODED  DATA* , 3X . afppop* ,4X , 

t/pcvo  data*,  ax,  /dfcoded  data*//) 

If  FORMAT ( SX. T4 ♦! ox, T1.13X, II ,11X, T 1 , 1 OX . 1 1 » 14X . T 1 ) 

17  FORMAT ( TSx • f 3) 

18  FORMAT ( l RX , 1 1 , 1 3X , T 1 , 1 1 X , 1 1 , 1 OX  » T 1 , 14* , 1 1 ) 

IQ  FORMAT ( 33* , J 1 , 1 1 X , T 1 , 1 OX . 1 1 ) 

20  FORMAT  (3X.*MonF=*»Tl  » 5X  , *S  YMD  ( * ",  T 1 , * > SI  =*  » II  , 2X  , *po=‘*  ,T  1 , 

1 2X  » *SB=* ,11) 

?1  FORMAT (5QX , 1 3 , *TW_R ANOOM  FRPQP  OQPPFCT IOM=* , T 1 ) 

’ 22  FORMAT (1X/*FRDOd"COUNT=XXX*/) 

23  FORMAT  <T3) 

24  FORMAT (1X/*PPP0R  POS T T I ON , F* , I 3 . * = x XX*/ ) 

25  format (Til 

c input  statfmrmts — 

PF  AD ( 2 » 1 7 ) LNOTM  

REAR (?, 17)  LR 
PFA0(R»17)  lx 

READ  ( 2 , 1 7 ) Tl.  AST  _ 

RF  AO  ( R • l 7 ) NR 
QFAD(?,17)  KR 
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RFA0(?»17)  LOCnO 
OFA0(?.17)  JPFH1 

RMP  (?»S)  JTM  

read  (?.5>  jpchr 

C TMTTI ALI7ATTOM  

LSYND=LX*LP»LMQTM  

ldlay=lp*lx 

MF=NR-KP 

JNpg  < TL  AST/YP ) _ 

LASTF= ( TL A5T/KP) *MP 
on  500  m=i,jmp 

DO  500  Nfl^NO 

j"RCVD(M.N)=0 
500  CONTINUE 

DO  510  K»1  »LMPTH  _ 

DO  510  L = 1 .kp 
KoEG(*.L>  =0 

5]Q  CONTINUE 

DO  SI?  L=1«KP 
Dn  51?  J=1.LDLAY 

SI?  JSTBOIJ«L>«0  _ _ 

00  514  N*T  .NS 
00  514  K=1 ,L5YN0 

p 1 4 J5  YND  ( K < U ) aQ 

no  5ia  t=i, Caste 

514  IFROR(1)=0 

WPITE (1.2?) 

RFAD  (4.??)  MM 
DO  3553  M=1*NM 

WP.TTF  ( 1 . 24  J M 

READ  (6.25)  T 
IEROR ( I ) =1 
3353  CONTINUE 
MODE=0 

00  520  M=1.MP 

DO  520  L=  1 .FP 

READ (?. 141 (KOFM (K,L»N) »Ksl .LNGTH) 

520  CONTINUE 

D0  530  JJ=l_*JPrHl 

DO  530  L = 1 

PFAO  (2.14)  (kORTM(k,|..N)  .K  = l .Lncoo) 

550  CONTTMUF 

WP  I TF  (3.0) 

on  540  l=i.nf 

DO  540  M=  1 . JPCH 1 

WP  I TF  (3.10)  m.L.  (F0RTH(k,L.u)  .k-  = 1 .L0O00) 
540  comttnuf 

C MOT«F  EMRFDDTMO  

CALL  onoFP 

1 = 1 

Do  6io  j=i  .jmp  ...  . 

DO  605  N=1 ,MP 


L 

I 
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JTfMPsNCOhr ( J,N) *TFO0O < f| 

JPCVD ( J.M) sMODTU ( JTFMP) 

1 = 1*1 

60*  COMTIMUF 
610  COMTIMUF 

c decoding 

DO  22??  Jsl.JNR 
DO  1500  M=l»NP 
IF  (M-KR)  1001*1001.1400 
1001  V=LX*L"P 
k=lngth 

JROUT  (N)=JSTRG(M.M) 

JOUT  ( J . M)  s.JPOI'T  ( SR 
1010  Ml=M-l 

JSTRG<m,M)=JSTRG(M1  ,M) 
r*“TMl)  1 0P0  * I 0P0  • 1 016 
1016  M=M-1 

GO  TO  1010 

~~l  ff?0— JS'TRG  («VM) VkOFG  (LNGTM.n) 

1030  K 1 =K- 1 

KPEG ( KjN ) =*RFG  ( K 1 , M ) 

IF  ( K 1 > ’1040.1040.100S 
1 00c  K=K-1 

GO  TO  1030  

10’40  KPEGff'.N)  a JPCVD ( J.M) 

GO  TO  1400 
1400  N1=N-KP 
JTFMP=0 
K=1 

L=  1 

bo  14?0  L=1.KP 
DO  1M0  Kal.LMGTH 

JTcMP=  KGFN(K.L.N)*KRFG(K.L)* ITFMP 

1410  COMTIMUF 
1 4?0  CONTTMtJF 

JTFMP=JTFmd* JOCVD ( J.N) 

JPGENINI ) =MOOTIJ  (.jTFMP) 

14Q0  GO  TO  1S0O 
1S00  COMTIMUF 

~C  PYMPPOMF-  P F G “oh  T F T I M6~TV N 0 ’ F S T I M A T I o N — 

Nl  = l 

DO  1640  Ml =1 .MS 
K=LSYMD 
1601  Kl =K- 1 

jsymd  (k,mh=joyno  <«■  1 .Ml  1 
IF  {(Cl  f 1 61 0*1610*1  GOT 
1607  K=K-1 

GO  TO  1601 

1610  JSYMD (K .Ml ) a JPGFM (Ml ) 

WR I TF  n.PO)  woor.Ml  .JSYMD  (1  .Ml  ) .JSYMD  (LDI  AY.Ml ) ♦ JSY-in  (lsymo.mI  ) 

1640__COMTTMUF 

IF  (MODF)  1POO,  1000.1000 


? 


* 


J 
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TMDFSHOI  D DECODFO  STMUL  STOP 


M.s.  * TM  PAGF  fi 


TiTOUB  r«02  L=1.KD 
Ll=L 

CALL  JROOM 

WRITF  <1.?1>  J.  JF(L) 

1«0?  CONTINUE 

EKLjaoo  L=1  <H'D 

IF  tjriu-l)  1890. 1820. 1«?0 
1820  DO  1880  M= 1 * MS 

Ml  =N»KO _ 

K1 =Ln*LNGTH-l 
00  1880  K = ?,|.NGTM 

JTEMPsJSYNO  (K1  «N)  ♦KGFM<K»N»M1  ) 

JGYND (K1»N) SMOOTH (JTFMP) 

K1=K1-1 

IPSO  CONTTMUF _____ 

1*80  CONTINUE 
1*90  CONTTMUF 

MODEs.Q _ 

60  TO  222? 

C R(J»ST  MOOF 

_JL90J2^  DO  1999  M=l,MF 

Ml =N+KR 

TF(JSYNOd.M))  1<302. 1OR0.1902 

1 °02  DO  1980  L=1.*P 

DO  19S0  Ks?«lmOTH 
Kl=LNOTH«.(_(a+LY-K*l 

IF_(KOEN(k.L,mJI)  ) 1911 .1012.101  1 

l9ii 

1012  GO  TO  10SO 

1 QSO  CONTTMUF  

JSYNO  i l ♦M>="0 
JF (L)=l 

JTEMPsJGTOG(LDLAY.L) *JF(L1 

j ’a  T P G ( L 0 L A Y . L ) = M O 0 T U { J T F m P ) 

1980  CONTTMUF 

1OO0  GO  TO  1990 

lOQQ  CONTTMUF 
2222  CONTTMUF 

C P o T N TOUT  QTatfmfmTS— - 

WPITF (3«IG> 

J=0 

JTN=J 

L = 1 

LAST  s JNP*MB 
Do  2100  T = 1 .1  AST 

!F1R1  ?O02, 2001 .2002 
2001  j=j*i 

WPITF  n.lS) J«  T NPUT  ( JTN)  « NOOOF  ( j«  1 ) ,TFOOP(T1  . 10CV0  ( 1*1  > .,I0HT|.I,I  > 

JTMsjfMVt 
GO  TO  ?0OQ 

2002_JF  0 -KP)  ?00C.?OOS,200S 

?npc  WPTTF  Cl. 18)  TNDUT(JTM)  .NCOOFf  I.L)  .TFOOPd)  « J0CVD(J.i  ) . JOi'Tf  J.L) 


r 


"i 
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JIN*J!NM 
GO  TO  ?0OQ 

2006  WPITF  (3. IP)  NCOOF ( J»L ) ♦ TFPOP ( T ) « JPCVD ( J»L ) 
IF  "(L-MP)  2007.20OR.200R 

2007  GO  TO  20PP 
200P  L=1 

GO  TO  ?100 
209P  L=L ♦ 1 

2J  00^  CONTI NUF  

CALL  FX IT 
FMO 
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APPENDIX  B 

DERIVATION  OF  ORTHOGONAL  EQUATIONS 


For  the  half-rate  convolutional  code  generated  by, 
G(D)  = 1 + D2  + D3  + D5  + D6  + D7  + D9  + D10  + D 


(1) 


parity  bits  are  generated  by  the  formula, 

P3+x+ll  ~ rB+x  + IB+x+l  + rB+x+2  + *B+x+4  + ^x+S  + rB+x+6 
+ 1 B+x+8  + IB+x+9  + XB+x+ll 

where  P denotes  a parity  bit,  I denotes  an  information  bit,  and  the  subscript 
notation  is  in  accordance  with  the  register  labels  used  in  Figure  Bl . 


(2) 


At  the  decoder,  shown  in  Figure  Bl  , a syndrome  is  generated  by  re- 
calculating parity  from  the  "noisy"  information  bits  and  comparing  this  parity 
bit  with  the  received  parity  bit.  The  current  syndrome  bit  is  thus  given  by, 

SB+x+ll  = ^B+x  + I’b+x+1+I'b+x+2  + rB+x+4  + rB+x+5 

+I  B+x+6  + 1 b+x+8  + 1 b+x+9  +I  B+x+11  + P B+x+11 

where  the  prime  (')  is  used  to  indicate  "noisy"  received  bits  which,  due  to 
errors,  may  not  equal  the  original  bits. 

Formally , 

rk  ■ \ * ** 


e\  • pk  * EPk 


(3) 


(4a) 
(4  b) 
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Figure  B-l . Biock  Diagram  - Adaptive 
Gallager  Decoder  (Kim's  Polynomial) 


Where  E_  is  0 or  1 according  to  whether  the  corresponding  received  bit 
is  in  error  or  not. 

Substituting  (4)  in  (3)  and  using  (2)  yields. 


SB+x+ll  ^B+x  + ^B+x+l  ^B+x+2  +E  B+x+4  +E  B+x+5 


+ E*. 


(5) 


B+x+6  +^B+x+8  B+x+9  +t  B+x+11  ^ B+x+11 


+ E 


+ E 


+ E1, 


From  Equation  5,  it  can  be  seen  that  the  values  of  the  syndromes  depend 
only  on  the  errors  in  the  information  and  parity  bits  from  which  they  are  formed 
and  not  on  the  actual  values  of  the  transmitted  bits  themselves.  Since  the 
syndrome  bit  is  formed  by  recalculating  parity  and  adding  i*  to  the  received 
version  of  itself,  it  is  clear  that  the  syndrome  will  be  zero  if  the  terms  in- 
volved are  all  correct. 

Conventional  threshold  decoding  makes  use  of  a set  of  orthogonal 
equations  generated  from  the  syndrome  sequence.  These  orthogonal  equations 
are  derived  from  linear  combinations  of  syndrome  equations.  These  are 
generated  in  the  following  manner.  First,  assume  that  there  were  no  transmission 
errors  prior  to  bit  position  X.  This  means  that  E^  =0  for  all  subscripts, 
j <x.  With  this  assumption,  the  following  set  of  syndrome  equations  can 
be  generated  from  Equation  5. 

s*  - ^ * eP* 

Sx+1  = e'x+I  + EPxtl 
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= El  + E1  . + E1  + EP 
x+1  x+2  x+4  x+4 

= Ei  +El  „+Ei  + E*  . + EP 

x x+2  x+3  x+5  x+5 

= E*  + E1  +E1  + E1  + E1  + EP 

x x+1  x+3  x+4  x+6  x+6 


E1  + E1  + E1  + E1  + E1  + El  + EP 
x x+1  x+2  x+4  x+5  x+7  x+7 

E1  . + E1  + E1  + Ei  + E*  E1  + Ep 
X+1  x+2  x+3  x+5  x+6  x+8  x+8 


Sx+9  = E1  + E*  + El  + E1  + E1  + E1  + E1  + EP 

x+Z  x+3  x+4  x+6  x+7  x+9  x+9 

Sx+10  = + £l  i + £l  „ + E1  + E1  + El  + E1  + E*  , . +EP 

x+iu  x x+1  x+3  x+4  x+5  x+7  x+8  x+10  x+10 

Sx+H  ■ ^x  + E‘x+l  + ^*2  + C4  * Cs  + * E‘X+8  *E‘x+9 

* E x+11  + ^x+ll 


The  following  set  of  orthogonal  equations  is  generated  by  forming  linear 

combinations  of  Equations  6.  The  subscripts  (with  the  x's  dropped)  indicate 

which  syndrome  equations  were  used  to  form  the  resulting  equation,  e.g. , 

A denotes  a linear  combination  of  syndrome  equations  Sx+4  and  Sx+7> 

**  / * 
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El  + EP 
x x 

E1  + E1  + EP  „ 
x x+2  x+2 


E*  + E1  + E1  + EP 
x x+1  x+3  x+3 


(7) 


M.7 


El  + EP  . + E1  + E1  +EP 
x x+4  x+5  x+7  x+7 


At  r Q = E*  + EP  . rP  . pi  + +E^ 

1'5'8  X x+1  + E x+5  E x+6  x+8  E x+8 


9,10,11  = Elx+Ex+4  +eFv+q  + E 


x x+4  x+9  x+10  + x+10  x+11 


+ e‘ 


x+11 


This  set  of  equations  has  the  orthogonal  properties:  (1)  Ei  appears  in 


every  equation,  and  (2)  no  other  bit  error  term  appears  more  than  once  in  the 


entire  set.  Since  the  set  is  orthogonal  on  E^,  it  is  possible  to  solve  for  this 
term  if  not  more  than  three  errors  are  present  in  the  terms  involved  in  the 
equations. 
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